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Целью работы является описание метода численного решения 
стационарного уравнения Шредингера, основанного на 
использовании интегрального уравнения, тождественного 
уравнению Шредингера. По своей структуре это интегральное 
уравнение близко к уравнению Фредгольма второго рода и 
допускает получение численного решения задачи. 
Рассматриваемый метод позволяет находить собственные 
значения энергии и собственные решения для квантово- 
механических задач различной размерности. Приведены 
результаты тестирования метода при решении задачи для 
одномерного и двумерного квантового осциллятора. Найденные 
численным методом значения собственных энергий и 
собственных функций осциллятора сравнивались с известными 
аналитическими решениями, затем оценивалась погрешность 
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результата. Наибольшая точность решения получена для первых 
энергетических уровней. Погрешность численного решения 
возрастает с номером собственного значения энергии. Для 
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последующего энергетического уровня погрешность 
увеличивается почти на порядок. Для четвертого 
энергетического уровня при разбиении области интегрирования 
на 500 элементов погрешность решения для собственных 
функций не превосходит 2%. Если энергетический уровень 
является вырожденным, то существует возможность получения 
всех собственных функций, соответствующих данному уровню. 
Для этого используются различные вспомогательные функции, 
симметрия которых согласуется с симметрией собственной 
функции. 
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фундаментальные решения. 


Введение. Фундаментальное значение при решении квантовых механических задач описывается уравнением 
Шредингера [1]. Аналитические решения этого уравнения могут быть получены лишь для весьма ограниченного круга 
задач, преимущественно одномерных. Поэтому разработано множество приближенных методов решения уравнения 
Шредингера, как аналитических, с использованием теории возмущения [2—4], так и прямых численных методов. 
Несмотря на широкий спектр имеющихся численных методов решения уравнения Шредингера, таких, например, как 
метод Нумерова [5], метод диагонализации [6-7], спектральный метод [7], и других численных методов [8-12], 
проблема эффективных способов нахождения собственных энергий и собственных функций для основного уравнения 
квантовой механики, особенно при решении многомерных задач, продолжает оставаться актуальной. Ниже 
предложена методика решения уравнения Шредингера, основанная на приведении его к интегральному уравнению с 
последующим численным решением. 
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Интегральное уравнение Шредингера. Пусть частица массой т, совершает финитное движение в силовом 


поле У (г). Запишем стационарное уравнение Шредингера в виде 
р УР р р 


Ат Еф = Ту). (1) 
й й 
Необходимо, решая уравнение (1), найти собственные значения энергии Ё и соответствующие им 
собственные функции, волновые функции частицы Ф () 
Перейдем от однородного уравнения Шредингера к неоднородному уравнению, введя новую искомую 
функцию 
О (г) =Ч(г)-щ (г), (2) 
где заданная вспомогательная функция ио(Г), как будет показано ниже, может определять некоторые особенности 
решения. Теперь уравнение Шредингера запишем как 


ДИ +0 =у(г)И — (А+А-+Ущ, (3) 
2т 2т 
где /, = 2 у(г) — РИ, 
Правую часть в (3) можно рассматривать как плотность некоторого фиктивного заряда р(г) ‚ создающего поле 


О (г) ‚ удовлетворяющее неоднородному уравнению Гельмгольца: 


д0+^ И =р(г), 
где 
р(г)=»(к)И (г)- (А+ +ь(г))ш (г). (4) 
Частное решение этого уравнения можно представить в виде 
р а (5) 


Здесь 2(т,К) — фундаментальное решение уравнения Гельмгольца, которое можно рассматривать как поле, созданное 
в точке г единичным положительным точечным зарядом, расположенным в точке К. Конкретный вид 
фундаментальных решений для уравнений Гельмгольца различных размерностей хорошо известен [13]. 

Физически разумно предположить, что область интегрирования (> в соотношении (5) близка, но несколько 
больше классически доступной для частицы области, границы которой определяются уравнением 

у (к) =: (6) 

В общем случае поле, соответствующее уравнению Гельмгольца, суммируется из поля, соответствующего частному 
решению (5), и поля, созданного внешними «зарядами», расположенными за пределами области @2, или на ее границе. 
Здесь также, исходя из физических соображений, можно полагать, что источники реального физического поля, поля 
волновой функции частицы, находятся в пределах области © и, таким образом, основная составляющая волновой 
функции И(г) удовлетворяет уравнению (5). Теперь, после подстановки в (5) плотности заряда (4), для искомой 
функции И(у) получается интегральное уравнение 


о а (®) (г, К)аО» (Т) 


Уравнение (7) есть интегральное уравнение Шредингера, являющееся следствием уравнения (3). Можно 
показать, что, наоборот, уравнение Шредингера (3) является следствием интегрального уравнения (7). 

Решение тестовой задачи для одномерного осциллятора. Уравнение (7) можно решать численно. Наиболее 
простой случай — решение одномерных задач. В этом случае, с помощью уравнения (6) находится отрезок [а, Б|= ©, в 
пределах которого потенциальная энергия частицы не превосходит ее полной энергии ЕЁ. На отрезке [а, 6] равномерно 
располагаются М узлов х; (1<1< М ). Значение поля в этих точках обозначим И;. Интегралы в уравнении (7) заменяем 
их приближенными значениями, используя какой-либо численный метод интегрирования, например, метод 
прямоугольников. В результате для № искомых величин И; получается система № линейных алгебраических уравнений 
вида: 

М 
>а0,=6. (8) 
1: 

При вычислении коэффициентов системы (8) при положительных /, используется фундаментальное решение 

одномерного уравнения Гельмгольца [14] 
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ех] вт (/%|х-Х|} 
в (хх) =—Щ о 
2./^. 
Описанный здесь метод численного решения уравнения Шредингера тестировался на примере решения 


задачи для квантового гармонического осциллятора. В этом случае область решения представляла собой отрезок 
ДЛИНОЙ 


2Е 
БЕ () 
то 


где « — частота осциллятора. Чтобы зафиксировать резкое убывание волновой функции за пределами классически 


ь 2Е 
допустимой области отрезок несколько увеличивают, полагая его, например, равным Ё=3,|—> . Затем задается 
то 


начальное значение энергии частицы, например, Р=0, и шаг варьирования энергии 0=. Желательно, чтобы этот шаг 
был много меньше расстояния ЛЕ между ближайшими собственными значениями энергии, 05<< АЕ. Для каждого 
значения энергии Е=ЁЕо+иб= последовательно, при п=0, 1, 2,... решается система (8) и находятся значения волновой 


функции в крайних точках области решения. Те значения энергии, для которых величинае = | |+ № м| имеет 


наименьшее. близкое к нулю значение, возможно, являются собственными значениями решаемой задач. Для 
уточнения этих значений можно СУЗИТЬ область поиска и уменьшить шаг варьирования энергии 05. Назовем параметр 


К =1/е=1/ [№] +|4 ^|) ‚ обратный к величине е, критерием задачи на собственные значения энергии. На рис. 1 


ь Е 1 ь 
представлен график зависимости критерия задачи К от приведенной энергии и=—--—, который в процессе 


йо 


вычислений менялся в пределах от 0 до 10 с шагом би = 0,1. Собственным значениям энергии в этом случае должны 


соответствовать максимальные значения критерия задачи К. Число узлов М в области решения задачи устанавливалось 
равным //=100. 
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Рис. 1. Зависимость критерия задачи на собственные значения энергии от приведенной 
энергии одномерного квантового осциллятора 


На рис. 1 наблюдаются четкие максимумы при целочисленных значениях и от 1 до 7. При и=0 максимум 
выражен достаточно слабо, а при и >7 положения максимумов несколько смещаются в область меньших значений и 
наблюдаются уже при и =7,9; 8,9; 9,8. Это, видимо, связано с тем, что с ростом номера собственной энергии 
количество осцилляций волновой функции в области решения возрастает. Это приводит к возрастанию погрешности 
численного решения для собственной функции и, соответственно, к увеличению погрешности для критерия К. Эта 
проблема устраняется путем увеличения количества узлов № и уменьшения шага варьирования энергии би. 


Численное решение уравнения (7) оказывается весьма чувствительно к значению параметра и, который 
определяет значение энергии Е. На рис. 2 представлены графики зависимости значений волновой функции \у от 
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приведенной координаты частицы лх/. Толстой сплошной кривой представлена собственная функция, 
соответствующая ц=2. График собственной функции, рассчитанный при количестве узлов М=100, в этом случае 
сливается с теоретической зависимостью, рассчитанной с использованием полиномов Эрмита. Пунктирная кривая на 
рис. 2 получена при и=2,1. 


в А 
РА 
г 





х/Г. 
Рис. 2. Зависимость значений волновой функции \у от приведенной координаты частицы х/ГЁ, 


Как видно из рис. 2, на графике, полученном при и=2.1, т. е. при значении энергии не равном собственному 
значению, кроме некоторого количественного различия с графиком собственной функции наблюдаются изломы, 
скачки. Наличие таких изломов, скачков является одним из признаков отклонения от собственного значения энергии. 

Основным параметром, определяющим точность численного решения, является, естественно, количество 
узлов № в области решения задачи. В качестве меры погрешности при вычислении собственных функций 
использовалась величина 





= = пах |; —$ф(х,) (10) 


где ф(х,) — значение аналитического решения задачи в точке с координатой х;. При вычислении погрешности 


численного решения предварительно производилась нормировка собственных функций путем их деления на 
максимальное по абсолютной величине значение этих функции в области решения задачи. На рис. 3 представлены 
графики зависимостей погрешностей численного решения от количества узлов М, полученные для трех собственных 
функций: с ц =0 (сплошная толстая кривая на рис. 3); и =1 (тонкая сплошная кривая на рис.3); и =2 (пунктирная 
кривая на рис. 3). 
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Рис. 3. Графики зависимостей погрешностей численного решения от количества узлов №. Сплошная толстая линия соответствует 
и =0; тонкая сплошная линия соответствует /=1; пунктирная линия соответствует д =2 
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Из рис. 3 очевидно, что погрешность численного решения для собственных функций убывает с ростом М, 
причем с увеличением № быстрота убывания погрешности уменьшается. Одновременно погрешность численного 
решения возрастает с номером собственного значения энергии: погрешность для последующего энергетического 
уровня увеличивается почти на порядок. Увеличивая число №, можно обеспечить приемлемую точность вычислений. 
Например, для четвертого энергетического уровня с и =4 при №М=500 погрешность решения не превосходит 2%. При 
этом следует отметить, что во всех случаях средняя или среднеквадратичная погрешность оказывается значительно 
меньше. 


Приведенные выше результаты получены при использовании вспомогательной функции и () =1+х. Если 


использовать вспомогательную функцию ис () =1, то при численном решении уравнения (7) правильные значения 
собственных энергий и собственных функций получаются только для четных значений и. Напротив, при 
использовании нечетных ц правильные значения собственных энергий и собственных функций получаются при 


нечетных и. Это обусловлено свойствами симметрии силового поля У(г)и фундаментальных решений, а также 


структурой уравнения (7). Используемая здесь функция из (Г) =1+х есть сумма четной и нечетной функций, что 


позволило получить решения, обладающие различной симметрией. 

Если при разработке алгоритма численного решения задачи учитывать определенные свойства симметрии 
волновых функций, можно оптимизировать компьютерную модель, сократить объем вычислений и повысить их 
ТОЧНОСТЬ. 

Решение тестовой задачи для двумерного осциллятора. Основным достоинством описываемого здесь 
метода решения уравнения Шредингера является возможность применения его, без существенного усложнения, при 
решении многомерных задач. Основная трудность решения связана с резким увеличением количества линейных 
уравнений типа (8) при увеличении размерности задачи. Ниже приводятся основные результаты решения тестовой 


двумерной задачи для двумерного изотропного квантового осциллятора. 
т 
При вычислении коэффициентов системы (8) при положительных значениях ^=—>Е используется 


фундаментальное решение двумерного уравнения Гельмгольца, выражаемое через функции Ханкеля нулевого 


порядка первого вида [13]: 
_ Ри = 
&(",К)= 4 Но (| К. 


Энергию осциллятора представим в виде Е=Ё,+Е,, где Е, и Е, — составляющие энергии, обусловленные 


движением вдоль осей Х или У соответственно. Область решения представляется в виде прямоугольника со 








28Р 2Е 
сторонами Д.. =3 - ир, =3 > . Разумеется, размеры области решения могут несколько отличаться от этих 
то то 


значений. Однако значительные отклонения этих размеров от рекомендуемых приводит к возрастанию погрешности 
вычислений, что может затруднить анализ результатов при определении собственных значений энергии. 

Алгоритм решения задачи на поиск собственных значений энергии повторяет описанный выше алгоритм для 
одномерного случая. Задается начальное значение энергии частицы, например, Ри=0, и шаг варьирования энергии 05. 


Для каждого значения энергии ЕЁ=Е\+п0& последовательно, при и=0, 1, 2,... решается система (8) и находятся значения 


волновой функции в крайних точках области решения. Те значения энергии, для которых величина е = Е м (здесь 
р. 


суммирование производится по значениям волновой функции | м в граничных узлах области решения) имеет 
наименьшее, близкое к нулю значение, возможно, являются собственными значениями решаемой задач. 
На рис. 4 представлен график зависимости критерия задачи на собственные значения энергии 


ы Е ы ы 
К =1/е=1/ + м от приведенной энергии и = =. —1, который в процессе вычислений менялся в пределах от —0,5 до 
Г © 


2,5 с шагом ой = 0,1. Собственным значениям энергии в этом случае должны соответствовать максимальные 
критерии задачи К. Число узлов М в области решения задачи устанавливалось равным М№ = 24х24 =576. Вычисления 


проводились при и () =1+х. 
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Рис. 4. Зависимость критерия задачи на собственные значения энергии от приведенной энергии 
двумерного квантового осциллятора 


На рис. 4 наблюдаются четкие максимумы при значениях и=0 и и=1. При и=2 максимум несколько смещен в 
область меньших значений, к и=1,9. Эти результаты согласуются с точными собственными значениями энергии 


двумерного осциллятора Ё,„ = 76 (и-+1), ци =0, 1, 2,... Однако следует обратить внимание также на наличие 


нескольких второстепенных локальных максимумов при значениях и, не совпадающих с теоретически ожидаемыми 
результатами. Увеличение числа узлов с №М=576 до М№=30х30=900 несколько улучшает результат, однако 
второстепенные максимумы сохраняются, что может приводить к ошибкам при идентификации собственных 
значений энергии. 

Определенные трудности возникают также при учете вырожденности энергетических уровней. Например, 
второму энергетическому уровню с и=1 можно поставить в соответствие две собственных функции. Одна из них 
симметрична относительно оси Х и антисимметрична относительно оси У, а другая — наоборот, симметрична 
относительно оси У и антисимметрична относительно оси Х. Для вычисления собственных функций в этих двух 
случаях описанным выше способом желательно использовать различные области решения (с разными соотношениями 


сторон) и различные вспомогательные функции из (г). В первом случае это может быть и (г) =1+ у, а во втором — 
Ид (г) Е 


Заключение. Проведенный анализ показывает перспективность предложенного численного метода решения 
квантово-механических задач различной размерности. Разработанный способ решения уравнения Шредингера может 
быть использован при решении различного типа эллиптических уравнений неканонической формы. Дополнительные 
возможности возникают также при использовании в интегральном уравнении, подобном уравнению (7), 
фундаментальных решений, предварительно найденных численным методом. Это позволяет, в частности, расширить 
круг решаемых задач путем использования метода точечных источников поля [15—20]. 
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